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Tins paper shows how C-numcrical-rangc related new strucures may arise from practical problems 
in quantum control — and vice versa, how an understanding of these structures helps to tackle hot 
vNtopics in quantum information. 

^ 1 We start out with an overview on the role of C-numerical ranges in current research problems in 
l/""^uantum theory: the quantum mechanical task of maximising the projection of a point on the unitary 
^f^prbit of an initial state onto a target state C relates to the C-numerical radius of A via maximising 
f^*3 the trace function | tr{C^U AU^}\. In quantum control of n qubits one may be interested (i) in having 
[I S SU(2 n ) for the entire dynamics, or (ii) in restricting the dynamics to local operations on each 
^■^qubit, i.e. to the n-fold tensor product SU(2) ® SU(2) (g • • ■ (g SU(2). Interestingly, the latter then 
J^^leads to a novel entity, the local C-numerical range W\ oc (C, A), whose intricate geometry is neither 
star-shaped nor simply connected in contrast to the conventional C-numerical range. This is shown 
^^ri the accompanying paper on Relative C-Numerical Ranges for Application in Quantum Control 
f-tind Quantum Information [1], We present novel applications of the C-numerical range in quantum 
'"ia;6ntrol assisted by gradient flows on the local unitary group: (1) they serve as powerful tools for 
"■yTlfeciding whether a quantum interaction can be inverted in time (in a sense generalising Hahn's 
f— f amous spin echo); (2) they allow for optimising witnesses of quantum entanglement. We conclude 
' I p y relating the relative C-numerical range to problems of constrained quantum optimisation, for 
C^wiiich we also give Lagrange- type gradient flow algorithms. 

• 'We are currently in the midst of a second quantum revolution. The 
, ^first one gave us new rules that govern physical reality. The second 



X 



one will take these rules and use them to develop new technologies. 

Dowling and Milburn, 2003 [2] 

_ i 

It was the main goal of the talk entailing this paper to entice the numeri- 
cal range community to showing further interest in problems of optimisation 
and control of quantum systems. By illustrating how important properties 
of C-numerical ranges relate to reachability and optimisation in quantum dy- 
namics, we wish to foster cross- fertilisation leading — in turn — also to new 
discoveries in numerical ranges. These may go beyond or follow earlier work 
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on applying higher-rank numerical ranges [3] to quantum error correction [4], 
or numerical ranges of derivations to anti-symmetric quantum states [5], on 
optimising coherence transfer in quantum systems [6-8], which entailed spe- 
cial interest in the C-numerical range of nilpotent matrices relevant in spec- 
troscopy [9, 10]. As contribution to our end, we present some new results on 
what we introduce as relative C-numerical range [1], we relate some of our 
recent results in quantum control to numerical ranges, and we express corner 
stones of quantum control in the setting of numerical ranges. A reader inter- 
ested in the quantum aspects may appreciate the paper being organised such 
as to pursue these issues in reverse order, whereas the one driven by impatient 
curiosity may prefer to jump into Section 11.51 or Chapters [2] and [3] right away. 



1 Overview on C-Numerical Ranges in Quantum Control 

Controlling quantum systems offers a great potential for performing computa- 
tional tasks or for simulating the behaviour of other quantum systems [11,12]. 
This is because the complexity of many problems [13] reduces upon going from 
classical to quantum hardware. It roots in Feynman's observation [11] that the 
resources required for simulating a quantum system on a classical computer 
increase exponentially with the system size. In turn, he concluded that using 
quantum hardware might therefore exponentially decrease the complexity of 
certain classical computation problems. Coherent superpositions of quantum 
states used as so-called 'qubits' can be viewed as a particularly powerful re- 
source of quantum parallelism unparalleled by any classical system. Important 
applications are meanwhile known in quantum computation, quantum search 
and quantum simulation: most prominently, there is the exponential speed- 
up by Shor's quantum algorithm of prime factorisation [14, 15], which relates 
to the general class of quantum algorithms [16, 17] solving hidden subgroup 
problems in an efficient way [18]. 

However, for exploiting the power of quantum systems, one has to steer 
them by classical controls such as voltage gates, radio-frequency pulses, or 
laser beams. It is highly desirable to do so in an optimal way, because the 
shapes of these controls critically determine the performance of the quantum 
system in terms of overlap of its actual final states with the desired target 
states. Here, the aim is to show how quantum optimal control relates to finding 
points on the unitary orbit of the initial quantum state A (in its density matrix 
representation) maximally projecting onto the desired target state C which is 
equivalent to maximising the target function 



F(U) = \tr{UArfCft}\ 



(1) 
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over all unitaries to yield the C -numerical radius of A. This is the scope in 
systems that are fully controllable in the sense that every propagator U £ 
SU(2 n ) can be realised on the physical system in question. 

Yet, often only local unitary operations K £ SU(2)®SU(2)®- ■ -<3SU(2) are 
actually available. Then the corresponding optimisation problems are confined 
to a subset of the conventional C-numerical range: this is the motivation 
to introduce the local C-numerical range W\ oc (C, A) as a special case of the 
relative C-numerical range Wk(C, A), whose intricate geometry is analysed in 
more detail in the accompanying paper [1]. 

In view of practical applications in quantum control, we finally give an 
outlook on contrained optimisation problems, i.e., those in which extremal 
points in the C-numerical range are searched subject to fulfilling contraints 
such as keeping UAU^ orthogonal to an undesired state D or leaving a neutral 
state E invariant. 

Since in general, a precise characterisation of solutions in algebraic terms is 
often beyond reach, we resort to numerical methods based on gradient flows 
on the unitary group. 



1.1 Quantum Dynamics: Notations and Relation to C -Numerical Range 

As usual in quantum mechanics, one may choose to represent a state of a 
pure quantum system by a state vector in Hilbert space, l^) £ "H- Its norm 
induced by the scalar product can be set to = 1, as will be assumed 

henceforth. The operators associated to quantum mechanical observables such 
as the Hamiltonian H are selfadjoint, so H = W. Then the Hamiltonian 
dynamics is governed by Schrodinger's equation of motion 

Mt)) = -%H m)) solved by \^{t)) = e"^ |^(0)> = U(t) |#))} . (2) 

Invoking Stone- von Neumann's theorem, the solution involves a time evolution 
by an element U(t) of the strongly continuous one-parameter group IA := 
{e~ ttH 1 t £ M} generated by the Hamiltonian H, see e.g. [19]. 

Non-pure quantum states comprise the settings of classically mixed states as 
well as reduced representations of quantum systems allowing for the descrip- 
tion of open dissipative systems. In these cases, one may choose to represent 
the state by a positive-semidefinite trace-class operator, the density operator 
p £ Bi(7i) with p > being normalised to tr p = 1. With the trace-class oper- 
ators Bi(H) forming a two-sided ideal in the bounded ones B(H), its dynamics 
follows Liouville-von Neumann's equation 



p(t) = -i[H,p(t)] solved by p(t) = U(t)p(0)U{tf , U(t) = e~ itH (3) 
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and thus 'dwells' on the unitary orbit of the initial state U (/9(O)). 

Next, consider the expectation value of observables B = G BiTL) in either 
setting. For pure quantum states it takes the form of the scalar product (•)•) 

(B)(t) := m)B\m) (4) 

while for non-pure states (using the Hilbert-Schmidt analogue tr{-' •}) it reads 

((B))(t) = tr (fit p(t) ) = tr ( B t u(t)p(p)U(tf) . (5) 

Clearly, in pure states the expectation value is an element of the field of val- 
ues (B)(t) G W(B) := {(u\Bu) \ \\u\\ = 1}, whereas in non-pure states it is an 
element of the C-numerical range {{B))(t) G W (B , p(0)) then taking the form 
of a real line segment. The latter is of particular significance, e.g., in ensemble 
spectroscopy, where it is customary to collect all signal-relevant components 
of the selfadjoint operator B in a new operator C that need no longer be 
Hermitian, and likewise the pertinent terms of p(0) in a general complex op- 
erator A. Thus moving from the selfadjoint operators B, p to arbitrary, not 
necessarily Hermitian (yet bounded) operators A, C, the analogue to the en- 
semble expectation value then becomes a general element of the C-numerical 
range W(C,A) := {tr(C t UAU^)\U G U{H)}, where here and henceforth 
A,C G Mat n (C ) are taken to be finite-dimensional in view of spin and pseudo- 
spin systems. 

The key features of W(C, A) [20,21] may thus be exploited for quantum 
optimisation and control. They comprise: (i) the C-spectrum of A is a subset 
of W(C, A); (ii) W(C, A) is always compact, connected and star-shaped [22] 
with respect to the centre tr{^4} tv{C^}/N; it is convex if (but not only if) C is 
normal with collinear eigenvalues in the complex plane, or if there is a p G C 
so that (C — pi) has rank 1; it is a circular disc in the complex plane [23] if 
there is a p G C so that (C — pi) is unitarily similar to block-shift form; (iii) 
the corners of the boundary dW(C, A) at which no tangent exists are always 
elements of the C-spectrum of A; (iv) for normal C with collinear eigenvalues 
in the complex plane (as well as in some degenerate cases [21]) W(C, A) is the 
closed convex hull of the C-spectrum of A thus forming a convex polygonal 
disc in the complex plane. 

1.2 Geometry of Optimisation within C- Numerical Ranges 

In the context of C-numerical ranges, there are geometric optimisation tasks 
immediately related to problems of quantum control, e.g., finding points on 
the unitary orbit of (the initial state) A that 
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(1) show a minimum Euclidean distance to (the target state) C corresponds 
to the maximum real part of the C-numerical range by 

\\C - U Atf\\l = \\A\\l + \\C\\ 2 2 - 2Retr{tf UAtf} j ( 6 ) 
while those that 

(2) enclose a minimal angle (mod ir) to C relate to the C-numerical radius 
r(C,A) := max|tr{C' t [A4[/ t }| by 

a | tr{Ctf7Af7t}|2 

cos AiC (U) = 2 2 . (7) 

Clearly, the mathematical limits to unitary transfer from A onto C are physi- 
cally meaningful only if all the transformations in the entire unitary group can 
be realised in the given experimental setting. This is what we will analyse in 
the next section. — In a fully controllable system, the C-numerical radius coin- 
cides with the maximal transfer of relevant components collected in A onto the 
target C. In coherent ensemble spectroscopy, this is identical to the maximal 
spectroscopic signal amplitudes obtainable in the absence of relaxation [6,7]. 



1.3 Controllability of Quantum Systems 

The standard bilinear control system with state X(t), drift A, controls Bj, and 
control amplitudes Uj G E reads 



X(t) = {A + Y j n J {t)B j )X(t) , (8) 
i=i 

while X(t) G QCn(C) and A,Bj G Mat j\r(C). Thus Hamiltonian quantum 
dynamics following Schrodinger's equation 

m 

|#)> = -i{H d + J2Mt)Hj) W)) (9) 

i=i 

m 

L>(t) =-i(H d + ^« i (t)fr i ) U(t) , (10) 

3=1 



can be cast into the above setting. Here ^ is the drift term, Hj are the control 
Hamiltonians with Uj(t) as control amplitudes. For n qubits, \yb) G C 2 ", 
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U G SU(2 n ), and %H V Gsu(2 n ). 
Definition 1.1 (Full Controllability) 

A system is fully controllable, if to every initial state A the entire unitary orbit 
O u (A) can be reached. 

In the special case of Hermitian operators this means any final state X(t) =: 
C can be reached from any initial state X(0) =: A as long as the operators A 
and C share the same spectrum of eigenvalues. 

Corollary 1.2 [24] 

The bilinear system of Eqn. Q3 is fully controllable if and only if drift and 
controls are a generating set o/su(2 n ) by way of the commutator, in the sense 
(H d ,H j \j = l,2,...,m) JjiB = Su(2 n ). 

Example 1.3 (n weakly coupled spin-^ qubits): 

Let a x = ( 5 J )> a y = ( ? ~o ) > ° z = ( o -l ) ^ e ^ ne P am i matrices. In n spins-|, 
a cifcz for spin fe is tacitly embedded as \ ®- -- \®a x ®1®---1 where a x is 
at position k. The same holds for ak y , ofe, and in the weak coupling terms 
°~kz°~tz with 1 < k < t < n. 

Theorem 1.4 (Controllability of Coupled Qubits) [7] 

A system of n qubits is fully controllable, if e.g. the control Hamiltonians Hj 
comprise the Pauli matrices {cr^x, o~k y \ k = 1, 2, ... n} on every single qubit se- 
lectively and the drift Hamiltonian encompasses the Ising pair interactions 
{Jki [pkzO~iz) /2 | k < £ = 2, . . . n}, where the coupling topology of Jk£ / may 
take the form of any connected graph. 

This theorem has meanwhile been generalised to other types of couplings 
[25,26]. 

Example 1.5 (Quantum Gates): 

In quantum computing, the logical gate operations have a unitary represen- 
tation. Therefore, implementing a unitary gate by a sequence of evolutions 
under drift and control terms of the respective hardware (i.e., the physical 
quantum system) can be seen as the quantum compilation task: it translates 
the unitary gates into the machine code of physically accessible controls. 

Corollary 1.6 The following are equivalent: 

(1) in a quantum system of n coupled spins-^, the drift and the controls 
Hj form a generating set of su(2 n ); 

(2) every unitary transformation U G SU(2 n ) can be realised by that system; 

(3) there is a set of universal quantum gates for the quantum system; 

(4) the quantum system is fully controllable; 

(5) the reachability set to the generalised expectation value ((C)) (t) = 
tr{C'A(t)} coincides with the C-numerical range W(C,A) for all A,C G 
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Mat 2 " (C). 

Proof : (1) 44> (2) follows by using the surjectivity of the exponential map 
in compact connected groups [24]; (2) 44> (3) express the same fact in the 
terminology of group theory (2) and quantum computing (3); (4) 44> (1) follows 
by Corollary 11.21 (4) 4^ (5) by definition of full controllability via reachability 
of the entire unitary orbit O u [A) . □ 



1.4 Tasks in Optimal Quantum Control 

'Optimise a scalar quality function subject to the equation of motion governing 
the dynamics of the system to be steered'-this is the formal setting of many 
an engineering problem both in classical and quantum systems. 

Extending the notions of quantum dynamics introduced above in Sec. 11.11 
there are two principle types of scenarios, (i) closed Hamiltonian systems 
evolving without relaxation and (ii) systems open to dissipation. Define 
U(t) := e~ ltH and the unitary conjugation map Adjj{-) : = C/(-)L r ~ 1 as well as 
the commutation operator adjj(-) := [H,-] to obtain the following equations 
of motion: 

(i) closed Hamiltonian systems 

1. pure state \j>) = -iH \i/>) \ip){t)eH = C N (11) 

2. gate U = -iH U U{t) G U{H) (12) 

3. non-pure state p = —i&dn (p) pit) G B\(Tt) (13) 

4. projective gate Adu = -i &d H ° Adu , Adt/(i) G U{Bi{H)) (14) 

(ii) open dissipative systems 

3'. non-pure state p = —(iadn +T) (p) p(t) G Bi(H) (15) 
4'. contractive map x = — (iad# +T) o x, x(t) G QL{B\(^H)j, (16) 

where Ti is a finite dimensional Hilbert space while U{7i) and B\(TL) denote 
the respective unitary group as well as the trace-class operators over TC, and 
QC(B\{TC)) is the general linear group over Bi{TL). Note that Eqn. [12] is the 
operator equation to Eqn. [11] referring to the unitary map of the entire basis. 
Likewise, for the unitary conjugation map, Eqn. [T3] is the operator equation 
to Eqn. [13] If the density operator p is viewed as a vector in Liouville space, 
e.g. by way of the vec representation [27], then the map Adu is an element of 
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the projective special unitary group 

(7(1) Zjv 

where Z at denotes the centre of SU(N). At the expense of being highly 
reducible, one may choose the embedded representation (U* <g) U) G SU(N 2 ) 
for the convenience of having (£/* (8) f7) vec(p) = UpU^ . 

When including dissipation via the positive semidefinite relaxation operator 
T > 0, the operator equation (Eqn. [TBI) to the Master equation (Eqn. [T5|) de- 
scribes the contractive quantum map generalising the unitary conjugation 
map Ad[/ for open dissipative systems. 

The scenarios of Eqns, IT2l 1141 1161 occur in the following typical optimisation 
problems of quantum control for (tracking over) fixed final times T: 

A Maximise Experimental Sensitivity in Coherent Spectroscopy by Finding 
the C -Numerical Radius of A [6,7]: 

maximise transfer amplitude f(U) := | tr{C"U(T)}| = | ti{0UAU^}\ 
subject to equation of motion U (t) = —iHU(t) 

B Realise Unitary Module Uq in Minimal Time: [28,29] 
maximise fidelity /(Adj/(T)) := Retr{Ad^ G Ad{/(T)} 
subject to equation of motion Adjj(t) = —i&dn ° Adjj(t) 

C Approximate Unitary Module Uq with Minimal Relaxative Loss: [30,31] 
maximise fidelity f(x(T)) := Retr{Ad^ x(T)} 
subject to Master equation of motion x(t) = — (i&dn +T) o x(t) 

Here we focus on problem A which determines the limit to unitary trans- 
fer on a general abstract level. It relates to the C-numerical radius in the 
fully controllable case and to the relative C-numerical radius [1] in the non- 
controllable case. In view of experimental implementation, in a second step, 
the family of critical unitary operators 

U :={U\\ tr(Ct UArf)\ = r c (A)} (18) 

may be realised or approximated in concrete experimental settings either in 
the fastest way (task B) or with least amount of relaxative loss (task C). 
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1.5 Gradient Flows Determining the C -Numerical Range and Radius 

In this section, we refer to numerical algorithms based on gradient-flows for 
obtaining the C-numerical radius, which means solving task A in the fully 
controllable case. 

If A, C G Mat n (C ) are Hermitian, the C-numerical range of A is a real line 
segment. Its maximum results from sorting the eigenvalues of A, C by magni- 
tude in same order as has been shown by von Neumann in 1937 [32] and in view 
of NMR spectroscopy by S0rensen [33]. For the special case of real symmetric 
matrices, a gradient flow on the group of special orthogonal matrices SO(N) 
was presented in the pioneering work of Brockett [34], a thorough analysis 
of which with convergence-ensuring step sizes for the numerical discretisation 
schemes can be found in the monography of Helmke and Moore [35] . While in 
the Hermitian case, the gradient flows generically converge to global extrema, 
an analogous result for the more general case, where A,C € Mat n (C) may 
be arbitrary complex matrices, is still missing, since it appears much more 
involved. However, the gradient flows may be generalised as shown in [6,8], 
and in all the cases we have been addressing over the years, the maxima found 
numerically have been on the boundary dW(C, A) as conjectured in [6]. 

With U := e~ itH and f(U) := tr{Ct UAW}, define the target functions 

F 1 ([/):= Re f(U) and F 2 (U) := \f(U)\ 2 . (19) 

For v = 1, 2 one finds the Frechet derivatives of F v at U £ U(N) as the linear 
maps on the tangent space TjjU(N) comprising elements of the form (iHU) 

DF u (iHU) = tv{G M iH} (20) 
G« = \{\UArt,Cft] - [UAU\C ] } ] ] (21) 

G< 2 ) = [u Art, eft] f(uy - [uArt,c^ f(u) (22) 

so gvadF u {U) = -G {u) U (23) 

is the respective gradient with G^ v ' skew-Hermitian. By compactness of U(N) 
the solutions of Eqn. [25] exist for all t € R and converge to the set of critical 
points since gradi 7 ^ is a real analytic gradient vector field [35]. 

Clearly, DF u (iHU) = in any direction H implies G^ v ' = 0. One may 
integrate the respective differential equation 



U = gvad F„(U) 



(24) 
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to arrive at the recursive scheme 

ui% = e~^U^ till ||G&|| 3 -0 • (25) 

As will be shown next, this gradient flow can readily be adapted to visualise 
the actual shape of W(C, A). 

Algorithm. (Determining the Boundary dW(C, A) [6,7]): 

The star-shapedness of the C-numerical range of A with respect to the centre 

tr A-tr /N [22] is central for the following straightforward gradient algorithm 

to determine the shape of W(C, A) by (best approximations to) its boundary 

dW(C,A): 

(1) shift the reference frame to the centre of the star: A i— ► A — tr ]y 1; 

(2) modify the above gradient algorithm to drive into the intersection of 
dW(C, A) with the positive real axis by the Lagrange approach described 
in the next paragraph; 

(3) rotate W(C, A) stepwise in the complex plane by way of multiplying say 
matrix A with a phase factor e li2w / m ^ 

(4) repeat steps (2) and (3) for I = 1, 2, . . . , m; 

(5) retransform the results in (2) into the original reference frame: the inter- 
section points then give an approximation of the circumference dW(C, A). 

Step (2) comprises a constrained optimisation implemented by Lagrange mul- 
tipliers as described next. 

Lagrange Approach to Tracing dW(C, A) [7]: 

Let again f{U) := ti{Cft UAU*} with U G U(N), where we assume without 
loss of generality the reference frame has been chosen such that A is traceless 
so the star centre coincides with the origin. In order to find the intersection of 
dW(C, A) with the positive real axis, one has to maximise F\(U) := Ref(U) 
while keeping F2(U) := Imf(U) zero. To this end, we introduced the Lagrange 
function 

L(U) := F X (U) - \{F 2 (U)f , (26) 
with A as multiplier. Its Frechet derivative has the components 



D{F!(U)} (iHU) = \ ti{{[UAU\C ] ] - [UAU\C^)iH} (27) 
D{F 2 (U)} 2 (iHU) = -i(F 2 {U)) ti{[UAU\&\ + [UAU\C^)iH} (28) 
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Figure 1. (Colour online) The C-numerical range for A,C of Example 1 1 . 71 using the Lagrange-type 
gradient-flow algorithm [7] described in the text. The dots give the C-spectrum of A, where the 
interior ones coincide with the star centre. 

so that one obtains the adapted recursive scheme [7] 

U k+1 = exp{-a {[U k AUlC ] ] s + 2i A (F 2 (U k ))[U k AUltf] H )} U k , (29) 

where for short, [•, -]s and [-,-]h denote the skew-Hermitian and the Hermitian 
part of the commutator, respectively. 

Example 1.7 [6] Define the following pair of 3 x 3 matrices 

A := diag ((0.7385 + 0.2400«), (0.0353 - 0.1660i), (0.4509 + 0.4060i)) and 
C := diag ((0.7189), (-0.1106 + 0.4789i), (-0.1106 - 04789i)) . 

(30) 

Fig. [T]shows the C-numerical range W(C, A) := tr{U AW C^} , where the con- 
cave triangle is a particular challenge to the algorithm, since it requires 500 
points to determine the circumference reaching into the vertices. The Lagrange 
parameter is dynamically increased with the iterations k [7]. Note that the 
vertex points derived from the C-spectrum of A are perfectly reached. Mean- 
while, the shape has been quantitatively corroborated by global optimisation 
methods — as has also been shown during the wonra in a collaboration with 
Prof. Tibken's group. 



2 Local C-Numerical Ranges 

In view of applications in quantum control, it is customary to term the n-fold 
tensor product SU(2) ® SU(2) (8) ■ ■ ■ <8> SU(2) as the group of local unitary 
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operations SU\ oc (2 n ). It is a subgroup to the full dynamic group SU(2 n ). 
Consequently, to a given initial state A, the reachability set under local con- 
trols amounts to the local unitary orbit 0\ 0C {A) := {KAK^\K £ SU\ oc (2 n )}. 

Definition 2.1 As in the companion paper [1], we define as local C -numerical 
range of A the subset 

W loc (C, A) := {tr (tfKAK*) | K £ SU loc (2 n )} C W(C, A). (31) 

It can be viewed as a projection of the local unitary orbit of the initial state 
A onto the target state C. 

As also seen in the concomitant study, in contrast to the usual C-numerical 
range, its local counterpart is no longer star-shaped, nor simply connected. 
With these stipulations, we will discuss recent applications of the local 
C-numerical range in quantum control. 

2.1 Application in Quantum Information 

Again, in terms of Euclidean geometry, maximising the real part in W\ oc (C, A) 
minimises the distance from C to the local unitary orbit 0\ 0C {A). 

In Quantum Information Theory, the minimal distance has an interesting 
interpretation in the following setting: let A be an arbitrary rank-1 state of 
the form A = \^){i>\ and let C = diag (1, 0, . . . , 0) £ Mat 2 »(C ). Thus in this 
case W\ oc (C, A) reduces to what we define as the local field of values W\ oc (A). 

Definition 2.2 (Pure-State Entanglement) 

An ra-qubit pure state A = \ip}(ip\ with \tp) £ C 2 " is termed a product state, if 
it can be written as a tensor product 

A = \lpl}(lpl\ <g> \i>2){^2\ <S> • • • <£> \lp n )(lpn\ with |^>£C 2 , (32) 

whereas it is said to be entangled if it cannot. 

Remark. In the present context, there are important observations with 
regard to the full unitary orbit O u (A) and the local unitary orbit Oi oc (^4) of 
pure states A of different types: 

1. all pure states form an equivalence class coinciding with O u {A) if A = 
IV'KV'I is an arbitrary pure state; 

2. generic elements on the full unitary orbit O u (A) of a pure product state 
A are pure, yet no longer of product form; 

3. all pure product states form an equivalence class coinciding with the 
local unitary orbit 0\ 0C (A) of an arbitrary pure product state A. 



Significance of C- Numerical Ranges and Local C -Numerical Ranges in Quantum Control 13 



Consequently, measures of entanglement remain invariant under local unitary 
transformation. 

Corollary 2.3 (Euclidean Measure of Pure-State Entanglement) 

For A = \tp)(tjj\, C = diag (1, 0, . . . , 0) £ Mat2^(C) the minimial Euclidean 

distance 

A:= min \\KAK^-C\\ 2 (33) 

KGSU loc {2") 

is a measure of pure- state entanglement because it quantifies how far A is from 
the equivalence class of pure product states. It relates to the maximum real part 
of the local numerical range W\ oc (A) via 

KAK* | \l = | \A\ \l + | |C| \\ - 2 Re tr{d KArf} 
= 2 — 2Retr{C^ KAK^} , 

where the last equality holds if also A is normalised to \ \A\\2 = 1. 

The (squared) Euclidean distance from the nearest pure product state is 
illustrated in Fig. [2] for the following two examples taken from quantum infor- 
mation theory [36,37]: 

Example 2.4 First, consider the pure 3-qubit state A(s) := | "03 ( s )) ("03 ( s )l 
parameterised by < s < 1 

\Ms)) ■= ^fs\W) + Vl^~s\W) , (35) 

where \W) := -^(0,1,1,0,1,0,0,0)* and \W) := 4g(0, 0, 0, 1,0, 1, 1, 0)* . 

Example 2.5 Likewise, observe the pure 4-qubit-state A(s) := \tp4(s)) {tp4(s)\ 



\Ms))-=V~s\GHZ , )-Vl^~s\^ + )®\^) , (36) 

where |GHZ 7 ) := -^(0,0,0, 1,0,0,0,0,0,0,0,0, 1,0,0,0)* and 
|^+)® \^+) := 1(0,0,0,0,0,1,1,0,0,1,1,0,0,0,0,0)* . 

Note that the plots in Fig. [2] quantitatively reproduce the global optimisation 
results found by numerical quadratic programming methods [36] yet cutting 
the CPU time by a factor of approx. 35 in Example 12.41 and by a factor of some 
150 in Example 12.51 Moreover, in both cases, the findings coincide with the 
exact solutions known algebraically [37]. 
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Figure 2. Euclidean distances of the pure states of Example 12.41 fa) and Example [23] (b) to the 
nearest pure product state as a function of s; here 'max. local transfer' denotes Retr{Ct KAK^}. 

Caveat: in non-pure states, the problem of entanglement is much more in- 
volved, since the state space to non-pure states forms no simplex: a generic 
density operator has infinitely many decompositions into pure states. For 
recent overviews, see, e.g., ref. [38,39]. Hence in those instances the above 
approach no longer applies. 



2.1.1 Significance of Entanglement. Finally, it is the purpose of this tuto- 
rial paragraph to show the numerical-range focussed reader why in entangled 
quantum systems, the total system comprises more information than accessible 
from putting together the information of all its local subsystems. 

To this end, we want to express a bipartite system (a, b) in terms of its 
subsystems a and b by making use of the respective orthonormal Hilbert space 
vectors \a u ) E TL a and G Tib. For simple demonstration, we assume both 
to be of the same dimension and consider a type of density operator in the 
total Hilbert space TL a <S> Tib that can be expanded in the Schmidt bases as 

Pab = 22 A At „(|a M )(a i y| <g> \b^){b v \) . (37) 

(IV 

The information locally accessible in subsystem a is encoded in the reduced 
density operator of subsystem a that is projected out by taking the so-called 
partial trace tr(|6 /1 )(6 !/ |) over the degrees of freedom of subsystem b yielding 

p a := tr b {pab} = tr 6 {E Ml/ A Ati ,(|a M )(a l/ | ® \bfj)(b v \)} 
= 5VW|a M )(a„| • tr{|6 M )(6„|}) 

(38) 

= S^A^(|a At )(a 1/ | • (Klbfj)) 
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where the last equality holds, because in the orthonormal Schmidt basis 
(bv\b/j,) = $fj,v The following standard examples [40] will illustrate reduced 
states in the scenarios of product states on one hand, and entangled states on 
the other. 

Example 2.6 Product states take the form p ab = p a <£> Pb, so one trivially 
finds 

p a = tr b {p a <8> Pb} = Pa • tr p b = p a (39) 

since tr p b = 1 by definition of a normalised density operator. 

Example 2.7 However, in the Bell state |$) := -^(|00) + |11)), where it is 

customary to use the short-hand |0) := (q) G C 2 and |1) := (5) as well as 
|00) := |0) ® |0) G C 4 and likewise |11) := |1) <g) |1), one obtains 

p a6 = |$}($| = i(|00) + |ll))((00| + (ll|) 

p a = tr 6 |$)($| 

= ±tr 6 (|00)<00| + |00)(11| + |11)(00| + |H)(H|) 

= i(|0)(0| o (0|0) 6 + |0)(l| a (l|0) 6 + |l)(0| a (0|l) 6 + |1)(1|„(1|1} 6 ) (40) 
= i(|0)(0| a + |l)(l| a ) 

j o 6 = tr a |$)($| = \t b 

The second example shows the generic situation: although p ab is a pure state, 
the reduced states of the respective subsystems p a ^ G T~£ a (b) are no longer 
pure. 

Remark. Compare the information content: 

(1) p ab = p\ h contains all information about the total system; 

(2) p a ( b j contains all information about the respective subsystem a(b); 

(3) the reconstruction p' ab := p a (8> p b ^ p'^ puts together all information 
accessible from both local subsystems; this is generically less than in p ab . 

In the second example, the reconstruction p' ab = p a p b = |l is diagonal, 
whereas the original p ab = Yl^u ^^v(\a^){a u \ ® \bfj,)(b u \) contained off-diagonal 
terms. Thus it is the coherent phase relation between the local constituents 
that is inevitably lost upon projection to the respective reduced systems. It 
cannot be reconstructed a posteriori using but local pieces of information. 
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This shows how in entangled quantum systems, the total system comprises 
more information than is accessible from putting together the information of 
all its reduced local subsystems. Measures of entanglement account for this 
loss of information and thus play an important role in quantum information 
theory. 

2.2 Application in Quantum Simulation 

In quantum control, it is of interest to decide, whether a given multi-particle 
quantum interaction (expressed by some interaction Hamiltonian H) can be 
sign-reversed solely by local unitary operations. If this is the case, one may 
undo or refocus the time evolution of such an interaction purely by local op- 
erations generalising the sense of the celebrated Hahn spin echo [41,42] in the 
following way: (0) start with any initial sate, (1) let the interaction evolve for 
some time t to give the propagator e~ ttH , (2) apply appropriate local opera- 
tions, (3) let the interaction evolve again for the same duration t, (4) apply the 
inverse to the previous local operations to (5) recover the initial state again as 
an echo, because steps (2)-(3)-(4) bring about the inverse propagator e +%tH . 
Note the same local operations apply to all the initial states; they only depend 
on the interaction Hamiltonian H. 
Mathematically, we ask whether there is a K € SU\ oc {2 n ) such that 

K e -itH R \ = e +itH fara ut G R+ . (41) 

Due to the equality K e~ ttH = e ~ ltKHK ' < the question readily boils down 
to deciding whether the sign-reversed Hamiltonian — H is on the local unitary 
orbit of the original Hamiltonian H, i.e., does there exist K £ SU\ oc (2 n ) such 
that 

K H K ] = -H . (42) 

Recently, we have solved this problem based on its normal form [43]: this is 
sign reversibility by local z-rotations, since every element K 6 SU\ oc (2 n ) is 
locally unitarily similar to local z-rotations. Here, we focus on the relation 
to local C-numerical ranges by illustrating that normalised Hamiltonians H 
are sign-reversible if and only if W\ oc (H, H) = [—1; +1], cf. [43]. To this end, 
define the spin-^ operators 

Jo:=l 2 ; Jz :=±diag(l,-l) ; J+:=(85); J- := 4 = (l o) ■ ( 43 ) 
In a single qubit, the J u G { Jq, J z , J+, J_} are the eigenoperators to the con- 
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jugation map Ad e -i#j z , i.e. 



Ade-^j, (J„) := eT^ 3 ' J u e 



Vi4>Jz 



-ip„<f> 



Jv 



(44) 



associated with the respective eigenvalues e ~ ip " ( ' } , po = p z = and p± = ±1. 
In order to generalise the arguments to the case of z-rotations on n qubits with 
individually differing rotation angles on each spin qubit <^>i, (fo, ■ ■ ■ , 4>e, ■ ■ ■ , 4>ni 
we write 



K 7 







-i<b?Jz 



-id n J z 



6 SU loc (2 n ). (45) 



Now consider a Hamiltonian in normal form H taking the special form 
H := H + + H_ with := S\, where H+ is a tensor product of Ad e -i«j z 
eigenoperators on each spin qubit I = 1, . . . , n according to 

H+ := JW jg) ® • • • <g> 4f ® • • • <g> (46) 

with independent z/^ € {0, z, +, — } on each spin qubit. With H± satisfying 



Ad 



(47) 



the Hamiltonian ^ is sign-reversed by local z-rotations provided there exists 
a set of rotation angles {4>{\ satisfying Yl!=iPvii ( t > i = i 71 " (mod27r) . This is 
the case if there is at least one spin qubit i giving rise to an interaction of 
quantum order p Ve = dbl. 



Moreover, a (real) linear combination 

m 



(48) 



A=l 



of Hamiltonians H\ := H\ + + H\_ with H\+ as in Eqn. 06] is jointly reversible 
by individual local z-rotations K z ((j>i, . . . ,cp n ) , if there is at least one con- 
sistent set of rotation angles {4>e} simultaneously fulfilling a standard linear 
system of m equations in n variables 



/ Piu! Piu 2 ■ 

P2V! P2u 2 ■ 
P3u t P3u 2 ■ 

\Pmv\ Pmv 2 



Plu„ \ / \ /7r(mod2vr)\ 

P2u n 
P3v n 



\4>n/ 



7T (mod 2-7t) 
7T (mod2-7r) 

\tt (mod 2tt) J 



(49) 
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With these stipulations, we have recently proven the following interrelations 
in view of local C-numerical ranges: 

Corollary 2.8 (Local Time Reversal [43]) 

For H = with \\H\\ 2 = 1 the following are equivalent: 

(1) the Hamiltonian H is sign-reversible under local unitary operations; 

(2) its local C-numerical range comprises —1, i.e., —1 E W\ oc {H, H); 

(3) its local C-numerical range is the interval [—1; +1] = W\ oc (H,H); 

(4) there exists a K 6 SU\ oc (2 n ) such that Adx(H) = —H; 

(5) H is locally unitarily similar to a H with Ad^ >( ^ n )(iJ ) = —H ; 

(6) H is locally unitarily similar to a linear combination of the form Eqn. 
satisfying the system of linear equations given in Eqn. Uty 

(7) let g = Qo © 0C£ , jj be the root-space decomposition of sl(N,C), where 

Eij denotes the square matrix differing from the zero matrix by just one 
element, the unity in the j th column of the i th row; then H is locally 
unitarily similar to a linear combination of root-space elements to non-zero 
roots (so i 7^ j) with Hy, '■= Xl™=i C A^' satisfying the system of linear 
equations Y17=l P-M ' fit = 7r(rnod 2tt) for X = 1, 2, . . . , m as in Eqn. \^9\ 

In ref. [43], we provided more tools to assess local reversibility by means of 
eigenspaces, graph representations of the interaction topology, spherical tensor 
methods, and root-space decomposition. Based on assertion (4) we also im- 
plemented a gradient-flow algorithm on the group of local unitaries SU\ oc {2 n ) 
in order to tackle the problem numerically. 

In the accompanying paper, we show the following: 

Theorem 2.9 (Local C-Numerical Ranges of Circular Disc Shape [1]) 
Let K be a compact connected subgroup of U(N) with Lie algebra t, and let 
t be a torus algebra oft. Then the relative C-numerical range Wn_{C,A + ) of 
a matrix A + G Mat tv(C) \ is a circular disc centered at the origin of the 
complex plane for all C € Mat at(C ) if and only if there exists a K G K and a 
A € t such that KA + K^ is an eigenoperator to adA with a non-zero eigenvalue 

&d A {KA + K r ) = [A, KA + K ] ] = ip KA + K ] and p + . (50) 

Clearly, ifKA + K^ is an eigenoperator of ad a to the eigenvalue +ip and A^ := 
A^ + , then KA_K^ shows the eigenvalue —ip. A + and A- share the same 
relative C-numerical range of circular symmetry, Wk(C, A + ) = Wk(C, A^). 

Moreover, locally reversible Hamiltonians H and nilpotent matrices 
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{A_|_,t4_} with rotationally symmetric local C-numerical ranges are related 
as follows: 

Corollary 2.10 

Let K = SU\ oc (2 n ). Let A + and A- := both share the same local 

C-numerical ranges W\ oc (C, A±) of circular disc shape for all C. 

Then 

(1) any linear combination A\ := A + + XA_ with A £ C and in particular the 
Hermitian H := A + + ^4_ are sign reversible by some local K £ SU\ oc (2 n ), 
so W\ oc (H,H) = [-\\A\\l;+\\A\\l], whereas 

(2) the converse does not necessarily hold: there are locally reversible Hermi- 
tian H to which no decomposition into a single pair {H + , H^} sharing the 
same rotationally symmetric local C-numerical range W\ oc (C,H±) exist, 
but 

(3) every Hermitian H £ Mat2™(C) that is locally sign reversible can be de- 
composed into at most ( 2 2 ) pairs (H+ , H^), (H+ , H^), . . . with each 
pair sharing the same rotationally symmetric local C-numerical range 
W loc (C,H { ^). 

Proof : 

(1) Eqn. [50] is equivalent to 

Ad e -*±(KA + rf) = e-^(KA + K^)e + ^ 

(51) 

= e"* adA (KA + K^) = e'^KA+K^ , 

thus p ensures a (f) with e ±J ^ p = — 1 to sign-reverse both KA±K^. 

(2) By Corollary 12.81 local sign-reversibility allows for linear combinations 
of eigenoperators to different ad^-operators with different eigenvalues 
thus generically violating the conditions for rotational symmetry of Theo- 
rem [2T9J 

(3) By Corollary 12.81 any Hermitian H £ Mat 2« (C ) locally reversible by 
z-rotations can trivially be decomposed into at most ( 2 2 ) Hermitian com- 
ponents (XijEij + \*-Eji) with l<i<j<N = 2 n and A^- 6 C , where 

the := (XijEij, X* j Eji) share the same rotationally symmetric 

local C-numerical range. □ 



3 Constrained Optimisation and Relative C-Numerical Ranges 



In quantum control, one may face the problem to maximise the unitary transfer 
from matrices from A to C subject to suppressing the transfer from A to D, 
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or subject to leaving another state E invariant. For tackling those types of 
problems, in ref. [7] we introduced a 'constrained C-numerical range of A\ 

Definition 3.1 

The constrained C-numerical range of A is defined by 



W(C,A)\ constia . nt := { tr(E/AE/t(7t) | constraint} C W(C,A) . (52) 



In ref. [7] we also asked which form it takes and — in view of numerical 
optimisation — whether it is a connected set with a well-defined boundary. 
Connectedness is central to any numerical optimisation approach, because 
otherwise one would have to rely on initial conditions in the connected com- 
ponent of the (global) optimum. 

Exploiting the findings on the relative C-numerical range of the accompa- 
nying paper [1], the structure of constrained C-numerical ranges can readily 
be characterised in some simple cases. 

Corollary 3.2 

The constrained C-numerical range of A is a connected set in the complex 
plane, if the constraint can be fulfilled by restricting the full unitary group 
U(N) to a compact and connected subgroup K C U(N). 

In this case, the constrained C-numerical range W(C, ^4) | constraint is identi- 
cal to the relative C-numerical range Wk(C, A) and hence the constrained 
optimisation problem is solved within it, e.g., by the corresponding relative 
C-numerical radius 7"k(C, A) . 

Proof : Direct consequence of the properties of the relative C-numerical range 
Wk(C, A) introduced in the accompanying paper [1]: if K is compact and 
connected, then Wk(C, A) is so as well, since it is a continuous image of a 
compact and connected set. □ 

Remark. Note that although being connected, Wk(C A) is in general neither 
star-shaped nor simply connected [1]. So if K is compact and connected this 
obviously extends to W(C, A) | constraint • 

3.0.1 Constraint by Invariance. The problem of maximising the transfer 
from A to C while leaving E invariant 




V 



(53) 



is straightforward in as much as the stabiliser group of E 



K E := {K £ U(N) | KEK^ = E} 



(54) 
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is easy to come by: it is generated by the Lie-algebra elements 

Ie '■= {k G u(N) | adjfc(-E') = [k,E] = 0} . (55) 

In particular, if E is of the form E = fil + Q with fj, G C and f2 G u (-^0) then 
fig is identical to the centraliser of in u(N). 

Lemma 3.3 

The set is closed under the Lie bracket, hence it is a subalgebra to u(N) 
thus generating a subgroup Kg C U(N), to wit the stabiliser group. 

Proof : Direct consequence of the Jacobi identity for the double commutator: 
[[A,B],E] + [[B,E],A] + [[E,A],B] = for all A,B,E G Matjv(C). Hence 
= [&2,-E] = implies A^], = and tg? is a Lie subalgebra to 
u(N) thus generating a compact connected stabiliser group K# C U(N). □ 

Remark. The stabiliser group of any E G MatTv(C) in U(N) is connected. 
This is in general not the case in SU(N) as easily seen for E := ({jj). However, 
one can restrict the above optimisation to the connected component of the 
identity matrix in SU(N) due to the invariance properties of the function 
U ■->■ tr{UAtftf}. 

A set of generators of t# may constructively be found via the kernel of the 
commutator map by solving a homogeneous linear system 

t E = kerad £ nsu(AO = {k G *u(jV)|(l<8> E- E 1 ® l)vec(fc) = 0} . (56) 

Corollary 3.4 

If solvable in an non-trivial way, the optimisation problem of Eqn. [53| entails a 
constrained C -numerical range that is connected. Moreover, it takes the form 
of a relative C -numerical range 

W(C,A)\ Adu(E)=E := {tr(UArfcft) | UErf = E} = W Ke (C,A) (57) 

and the optimisation problem is solved by the relative C -numerical radius 
rK B (C, A). In Hermitian E, ~Ke includes a maximal torus group T C SU(N). 

Proof : To any E G Mat7v(C ) the existence of a compact connected stabiliser 
group C U(N) can constructively be checked as in Eqn. [56l Moreover, 
every Hermitian E can be chosen diagonal: hence in that case Eg includes a 
maximal torus algebra t with tC ff; C u(N). The rest follows. □ 

Remark. Obviously, the constraint of leaving E invariant while maximising 
the transfer from A to C only makes sense, if A and E do not share the same 
stabiliser group. 
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Since it may be tedious to check for the stabiliser group of E in each and 
every practical instance and then project the gradients onto the corresponding 
subalgebra %e, a more versatile programming tool would be welcome. 

Algorithm. To this type of constrained optimisation, in ref. [7], we also derived 
a gradient flow based on the Lagrange function (with fc(U) '■= tr{C^C/ 'AU^}) 

L(U) = \f c (U)\ 2 - A (tr{UErtrf} - , (58) 

where the constraint U EU^ — E = was written in the more convenient form 
tx{UEU*E* } — \\E\\\ = 0. The algorithm implements the gradient from the 
Frechet derivatives 

D{\f c (U)\ 2 - Xf E (U) + A (iHU) 

(59) 

= tr { (2 (fc(U)[UArf,Cft]) s - \[UEU\E^}) iH) , 
where for short, (-)s denotes the skew-Hermitian part, within the recursion [7] 
U k+1 = exp { - a (2 {fc(U k )[U k AUlc% - \[U k EUl^\) } U k . (60) 



3.0.2 Constraint by Orthogonality. For the sequel of this paragraph, we 
assume a generic shape of the C-numerical range W(D, A) so that uiq G 
W(D, A) defines the unique point in W(D, A) that is closest to the origin in 
the complex plane. For these instances, we address the optimisation problem 

max\ti{UAU j C j }\ subject to tr{ U AU ] D ] } = m . (61) 

Clearly, perfect matches exist if £ W(D,A), because only then are there 
points on the unitary orbit O u (A) that are orthogonal to D. Moreover by 
Eqn. [71 the modulus of tjiq relates to the cosine of the angle between D and 
points on the unitary orbit O u {A) that come closest to orthogonality. The 
corresponding C-numerical range constrained by (best approximation to) or- 
thogonality to D takes the form [7] 

W(C,A)\ Ad ±D := {tiiUArfCft} | \ tx{UAU^D^}\ = m } C W(C,A) , 

(62) 

which is more difficult to characterise, because in order to establish whether 
it is a connected set in the complex plane, one has to check the constraint set 



K D := {K C SU(N) \ tr{KAK^D^} = m } 



(63) 
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for its properties of (i) forming a subgroup and (ii) connectedness. Gener- 
ically, already the first condition is violated, and only in the rare event of 
both of them being fulfilled, the set Kd turns into a subgroup Kd, and 
hence the constrained C-numerical range into the relative C-numerical range 
Wk_ d (C, A). Then the constrained optimisation problem of Eqn. [B"T1 would 
be solved by the relative C-numerical radius /"k d (C, A). Even in the special 
case rriQ = £ W{D, A) it is difficult to make sure the projection of the 
unitary orbit O u (A) onto the orthocomplement D 1 - of D in the Hilbert space 
(Mat tv(C ), tr{-'-}) is still a smooth manifold. Generically, again this is not 



It is for these reasons that addressing orthogonality problems by a Lagrange 
approach is more promising. For this to make sense, one trivially has to ensure 
C and D are not scalar multiples of one another, yet is not necessary that 
perfect orthogonality in the sense of 6 W(D, A) can actually be achieved. 

Algorithm. In ref. [7], we devised a Lagrange- type gradient-flow algorithm for 
solving the constrained optimisation problem of Eqn. IBTI numerically. To this 
end, define f c (U) := tv{UAU^C^} and f D (U) := tt{UAU^D^}. Introducing 
the Lagrange function 



the task to maximise the transfer from A to C while suppressing the transfer 
from A to D can be addressed by implementing the gradient of 



D{\f c (U)\ 2 - X\f D (U)\ 2 } (zHU) = 

tr{2(ft(U)[UArt,Cft])sin} - A tr {2 {f D {U)[UAU\ D*]) s iH} 



U k+1 = exp{-2a{(r c (U k )[U k AUlC j }) s - A (fh(U k )[U k AUl &]) s )}U k 



the case. 



L{U) :=\f c (U)\ 2 -\\f D {U)\ 2 



(64) 



into the recursive scheme [7] 



(66) 



Example 3.5 In Fig. [3l a first example is given for the matrices 




0.8359-0.1152i 







-0.2593-0.3906* 







0.0151+0.2609i 



) 



(67) 




-0.0318+0.0690i -0.3522-0.3185i 0.2351-0.3050i 
-0.0404+0.0656i 0.0844-0.2880* 0.2135+0.3234i 
0.3086+0.1076* 0.1742-0.2291* -0.2368+0.3585* 



) 



(68) 



/ -0.2910-0.3480* -0.2395+0.0274* -0.2428+0.0656* \ 
D = 0.0836-0.2790* -0.1836-0.0203* -0.2427+0.2396* 
V -0.3906-0.1387* 0.1989-0.2725* -0.0442+0.3871*/ 



(69) 
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(a) (b) 




-1 D 1 -1 1 

Re Re 

Figure 3. Constrained optimisation A — > C suppressing A — > D for matrices A, C, D of Ex. 13,51 

the first one of which is normal thus entailing a triangular pattern in its nu- 
merical range VF(^4) shown together with W(C') and W(D^) in column (a). 
The conventional numerical ranges were calculated using the classical algo- 
rithm of Marcus [44-46]. By the constraint that the transfer from A to D 
be minimal, the local maxima on the boundary dW(C, A) of the C-numerical 
range of A (asterisks in column (b)) are shifted to the points indicated by 
noughts. Note that only the upper left one seems to be displaced from the 
boundary dW(C, A) slightly into the interior of W(C,A). All the 100 trajec- 
tories (shown as dotted lines) starting from random initial conditions converge 
into the same maxima, while the transfer from A to D gets zero as shown in 
the C-numerical range W(D, A) at the bottom of column (b). 



4 Conclusions 

We have shown how C-numerical ranges provide the setting for many quantum 
optimisations. Knowing about its structure paves the way to numerical algo- 
rithms, e.g., for plotting its shape. In the accompanying paper, we introduced 
the relative C-numerical range of A as a restriction of the full unitary group U 
to some compact connected subgroup K C U, in which case it is connected but 
not simply connected. Here we illustrated that this is of practical importance 
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in relevant examples from quantum information: For instance, the maximum 
real part of the local C-numerical range of A (as a special case of the relative 
C-numerical range) directly corresponds to a measure of 'pure-state entan- 
glement'. In other instances, if the local C-numerical range of a normalised 
interaction Hamiltonian H is W\ oc (H, H) = [— 1,+1], then the interaction is 
locally reversible. These cases are fully characterised in Lie algebraic terms 
(by root-space elements to non-zero roots fulfilling a linear system of equa- 
tions) and they are related to cases, where the local C-numerical range is a 
circular disc in the complex plane. 

Moreover, some constrained optimisation problems in quantum control can 
be treated group theoretically: if the constraining conditions can be trans- 
lated into restricting the full quantum dynamics on [/ to a compact connected 
subgroup K, then the optimisation problem remains within the corresponding 
connected relative C-numerical range and the optimisation amounts to finding 
its relative C-numerical radius. For more general cases we provide numerical 
algorithms for constrained optimisations of Lagrange-type. 

5 Outlook 

Motivated by applications in quantum control, the relative C-numerical range 
of A introduced awaits further mathematical elucidation: for instance, what 
are the properties of its boundary, under which conditions is it simply con- 
nected, or even star-shaped? When is it convex? How can one systematically 
find compact connected subgroups embracing practical constraints of quantum 
optimisation so that they relate to a connected relative C-numerical range? 
What happens in generalisations where the subgroups are no longer compact 
and connected? Are there simple instances, in which the correponding re- 
stricted C-numerical ranges have few connected components and thus are not 
'hopeless' in view of practical optimisation? Can one prove the conjecture of 
ref. [6] that in generic cases, the gradient flows of Section [1.51 always converge 
to points on the boundary dW(C, A) and there are no local maxima in the 
interior of W(C, A)? 

We anticipate that problems in quantum control will profit from mathematical 
results addressing those questions, and — in turn — studying quantum dynam- 
ics will stimulate conceiving new structures worthy of mathematical research. 
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